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Nuclear magnetic resonance (NMR) has been widely used in the graphene oxide (GO) structure 
study. However, the detailed relationship between its spectroscopic features and the GO structural 
configuration has remained elusive. Based on first principles ^^C chemical shift calculations using 
the gauge including projector augmented waves (GIPAW) method, we provide a spectrum-structure 
connection. Chemical shift of carbon is found to be very sensitive to atomic environment, even with 
an identical oxidation group. Factors determining the chemical shifts for epoxy and hydroxy groups 
have been discussed. GO structures previously reported in the literature have been checked from 
the NMR point of view. The energetically favorable hydroxy chain structure is not expected to be 
widely existed in real GO samples according to our NMR simulations. The epoxy pair we proposed 
previously is also supported by chemical shift calculations. 



I. INTRODUCTION 

Graphene has been the rising-star material in the past 
several years, due to its unique electron structure and 
its great application potential. 1-4] However, its large- 
scale production remains a big challenge. Reduction of 
graphene oxide (GO), which is obtained by oxidation and 
exfoliation of graphite, Q provides a promising way to 
obtain graphene or chemically modified graphene mas- 
sively. |6|-l9| Solution based chemical processes are ex- 
pected to be applicable to many kinds of applications, 
which makes GO itself an important material. (Tol - [l5| 

Despite a big effort from both experimentalists [1, [l6l - 
[2l| and theoreticians, [23 - [27| the precise chemical struc- 
ture of GO is still not clear. The main reason for the 
difficulty to understand its structure is the amorphous 
and nonstoichiometric nature of GO. At the same time, 
different samples have different levels of oxidation, which 
makes things even more complicated. Nevertheless, some 
fundamental structural features of GO, as proposed by 
Lerf and coworkers, [l^ have been widely accepted. In 
this so-called Lerf mode, hydroxy (-0H) and epoxy (- 
0-) groups spread across the graphene planes, while car- 
boxylic acid groups (-COOH) exist at edge sites, possibly 
in addition to the keto groups. 

Aside from this general picture, there are still many 
open questions about the GO structure. For example, 
how are the hydroxy and epoxy groups distributed on 
the GO plane? Do they aggregate together or avoid each 
other? Is there any other new species in GO, especially 
for highly oxidized samples? Are sp^ carbon atoms clus- 
tered in aromatic forms? 

It is very difficult to answer these questions based only 
on experimental raw data. Therefore, theoretical model- 
ing has played an important role in GO structure study. 
Based on first principles energetics, many GO structure 
models have been proposed. Adapting a 2x2 unit cell. 
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Boukhvalov et al. |2J] and Lahaye et al. [l^] system- 
atically studied the possible atomic configurations with 
different ratios of epoxides to alcohols at different degrees 
of oxidation. Based on a study with similar strategy, Yan 
et al. (25| emphasized that it is favorable in energy to 
form hydroxyl chains with hydrogen bonds. In a GO hy- 
drogen storage study, Wang et al. ,121 constructed a GO 
model with a building block which has the lowest energy 
at a very low oxidation group coverage. 

Although it is very straightforward to conduct in theo- 
retical study, the power of energetics analysis is expected 
to be limited by the complexity of the GO potential en- 
ergy surface, especially when artificial periodic boundary 
condition must be adopted with a small unit cell. In 
contrast, computational spectroscopy provides informa- 
tion which can be directly compared with experiments. 
Therefore, it is a powerful alternative in nanostructure 
studies. For example, using Raman spectrum simula- 
tion, Kudin ct al. [28] proposed an alternating single- 
double bond GO structure model. Recently, our x-ray 
photoelectron spectroscopy (XPS) simulation also led to 
a prediction of the epoxy pair and epoxy- hydroxy pair 
species in highly oxidized GO samples. |26| 

Nuclear magnetic resonance (NMR) is the most widely 
used experimental technique in GO structure study. 
[9l[l6l-[20j In fact, the popular Lerf model is mainly based 
on ^"^C NMR experiments, [l^ Currently, very high qual- 
ity NMR spectrum is available via synthesis of almost 
fully i^C-labeled GO. [13] There are mainly three broad 
resonances in the ^'^C NMR spectrum of GO. By check- 
ing a series of reactions of GO with different reagents, the 
peak around 60 ppm is assigned to carbon atoms bond- 
ing to the epoxy group, and the peak around 70 ppm 
is corresponding to the hydroxy group connected carbon 
atoms. [13] The sp^ carbon has a NMR peak at about 
130 ppm. Besides the three main peaks, there are also 
three small peaks clearly shown in the high resolution ^^C 
NMR spectrum. [2^ Gao et al. [Qj suggested the these 
three peaks at 101, 167, and 191 ppm are corresponding 
to lactol, carboxy, and ketone groups, respectively. All 
these NMR assignments are based on chemical intuition, 
it is highly desirable to confirm them by first-principles 
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simulations. 

Theoretical NMR simulation has been widely used in 
studies of carbon nanotubes (CNTs). (H-HI ^^C NMR 
isotropic chemical shifts of pristine CNTs, based on the 
density functional theory (DFT) and an infinite CNT 
model, depend on the electronic structure and the di- 
ameter but not on the chirality. [sH] It has also been 
demonstrated that a wealth of knowledge on functional- 
izations [H,!!!] and defects [111 of CNTs can be obtained 
from NMR simulations. Similarly, new insights about 
GO structure are expected to be obtain by a systematic 
computational NMR study. 

II. COMPUTATIONAL DETAILS 

All of our calculations were performed with a plane 
wave based DFT implementation within the CASTEP 
package, [H, [S^ and the gauge including projector- 
augmented plane-wave (GIPAW) method [33, |40| was 
used to calculate NMR shielding tensors. Ultrasoft pseu- 
dopotentials [13] were used to describe the interaction 
between valence electrons and ions. The revised Perdew- 
Burke-Ernzerhof (RPBE) density functional [s^ used in 
this study has been demonstrated to be accurate enough 
for NMR simulation of CNTs. [3(| The energy cutoff of 
plane wave basis set in our calculations was chosen to be 
high enough (500 eV) to converge energy, geometry, and 
nuclear shielding tensor. 

When considering isolated oxidation groups, to avoid 
handling metallic system, we adopted an armchair 
graphene nanoribbon (AGNR) with a finite band gap as 
our substrate system. The vacuum space between two 
neighboring AGNRs was wider than 10 A in both direc- 
tions perpendicular to the ribbon direction (the z direc- 
tion in Figure [1^) . Along z direction, unless otherwise 
specified, a two-unit supercell (8.52 A long according to 
the graphene lattice constant) was used to make sure 
that oxidation groups on the ribbon do not interact with 
their images in neighboring supercells. Atom coordinates 
of the functionalized ribbons were fully optimized while 
the supercell size was fixed. A 1 x 1 x m fc-grid was used 
for Brillouin zone integration, and m was large enough 
to make the isotropic chemical shifts of all carbon atoms 
differ less than 1 ppm to those with a lxlxm-1 grid. 




FIG. 1: An epoxy group adsorbs at (a) a P-site or (b) a D-site 
on 9- AGNR. (c) and (d) Two full-oxidized GO with epoxide 
only. Supercells used in our calculations are marked by red 
line. Carbon is in gray, oxygen in red, and hydrogen in white. 



For highly oxidized GO samples, graphene based two 
dimensional (2D) models were used. And some special 
structure models proposed in previous energetics studies 
were taken for NMR calculations. We optimized both 
atom coordinates and cell parameters for all 2D models. 
The fc-grid was also carefully tested to make sure that the 
NMR chemical shifts of all carbon atoms have converged. 

As a commonly used method, the calculated ^'^C chem- 
ical shift of benzene was used to calibrate the chemical 
shift for other systems, [l^llO] which generally leads to a 
quantitative agreement between theory and experiment. 
The isotropic chemical shift of benzene molecule calcu- 
lated using a cubic cell of 15x15x10 A'^ was 42.2 ppm, 
agreeing well with previous studies, [s^j As an additional 
test, we calculated the chemical shift of a (7,0) CNT, and 
the result, 136.3 ppm, also agrees well with the value 
(136.4 ppm) obtained by Zurek et ah. [s^l 



III. RESULTS AND DISCUSSION 
A. Epoxy groups 

First, we consider an isolated epoxy group on an 
AGNR with 9 carbon chains (9- AGNR) . The epoxy group 
can be attached to a C-C bond either parallel to the rib- 
bon direction (P-site, Figure [T^) or diagonal to the rib- 
bon direction (D-site, Figure [T}d). Cb is used to name 
carbon atoms connecting to an epoxy group. The op- 
timized C-C bond length between the two C^; carbon 
atoms is 1.50 A for the P-site and 1.52 A for the D-site, 
close to that of sp^ hybridized carbon. For a P-site ad- 
sorbed epoxy group, chemical shifts of the corresponding 
two Ce atoms are both 68.0 ppm. [4l| Due to the its 
lower symmetry and the finite AGNR width, there is a 
small difference (less than 1.5 ppm) between the chemi- 
cal shifts of the two Ce atoms in the D-site case. Their 
averages is 69.4 ppm. It is very interesting to note that 
the chemical shift of isolated epoxide group is close to 
the experimental peak value of hydroxy group around 70 
ppm instead of 60 ppm assigned for epoxy group. 

To check the width effect of the ribbon model, we 
have also tested wider AGNR with 10 carbon chains (10- 
AGNR). The corresponding C-C bond lengths are 1.50 
and 1.53 A for the P- and D-site, respectively. The av- 
erage chemical shifts of the two C^ atoms are 67.2 and 
73.2 ppm for these two adsorption sites. Their difference 
becomes larger than that of 9-AGNR, which may due 
to the 3-family behavior of the ANGR electronic struc- 
ture. [13 However, the overall picture for 10-AGNR is 
the same as before, and chemical shifts close to 70 ppm 
are still obtained for individual epoxy groups. 

It is interesting to see how the chemical shift changes 
upon the change of the epoxy group concentration. As 
another limit compared to the isolated epoxide case with 
the lowest concentration, we consider the fully epoxidized 
GO, which has the highest epoxy group concentration. 
Two fully-oxidized models with chemical formula C2O 
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have been considered. Both have an optimized C-C bond 
length of about 1.50 A. The first one (Figure[TJ;) has a 
chemical shift of 59.8 ppm. The second one (Figure [Iji) 
is part of a previously proposed GO structure model, [25] 
which leads to a Cb chemical shift of about 57.4 ppm. 
Therefore, the 100% functionalized graphene epoxide has 
a chemical shift close to the experimental value 60 ppm, 
which is much lower than that of isolated epoxides. ^'^C 
chemical shift of epoxide in GO is thus very sensitive to 
chemical environment, and isolated epoxy groups should 
not be widely existed in the experimentally prepared GO 
samples. 

B. Hydroxy groups 

For single hydroxy groups, a four unit supercell along 
the ribbon direction with two hydroxy groups has been 
used to calculate their chemical shifts. j41|| Ch is used 
to indicate carbon atoms connected to a hydroxy group. 
For 9-AGNR (Figure [2}i) , the chemical shift of the 
isolated hydroxy groups is 72.4 ppm, close to the 70 ppm 
peak position in experiments. 10-AGNR gives a similar 
result, with a Ch chemical shift of 72.2 ppm. 

Isolated hydroxy pairs have also been considered. En- 
ergetically, the most stable configuration has two hy- 
droxy groups occupying a 1,2-site on two opposite sides 
of the graphene plane. P-site and D-site can be simi- 
larly defined as in the isolated epoxy group case. For 
9-AGNR, the chemical shifts of the two Ch atoms are 
very close, and their averages are 72.2 and 72.0 ppm in 
the P-site (Figure [2)d) and the D-site (Figure [2};) cases, 
respectively. Both values are close to 70 ppm assigned to 
hydroxy groups in experiment. The optimized Ch-Ch 
bond length is around 1.54 A for both the P- and D- 
sites. Similar results are obtained for wider 10-AGNR, 
and the averaged Ch chemical shifts are 68.2 and 73.3 
ppm for the P-site and D-site, respectively. 



FIG. 2: (a) An isolated hydroxy group, (b) a P-site hydroxyl 
pair, (c) a D-site hydroxyl pair, (d) a 1,3- hydroxyl pair, (e) 
a 1,4-hydroxyl pair on the same side, and (f) a 1,4-hydroxyl 
pair on two opposite sides. Only part of the supercell used 
in calculations is shown. [4lj| Carbon is in gray, oxygen in red 
and hydrogen in white. 

To further understand the chemical shift of isolated 
hydroxy pair, three other models have also been consid- 



ered, where the two hydroxy groups occupy a 1,3-site on 
two opposite sides (Figure [2ji), a 1,4-site on the same 
side (Figure [2I3) , or a 1,4-site on two opposite sides (Fig- 
ure Hf). The 1,3-hydroxyl model with hydroxy groups 
on the same side is not stable after optimization. For 

9- AGNR, the chemical shifts of the two Ch atoms in the 
opposite-side 1,3-hydroxyl case is close, and their aver- 
age is 70.8 ppm. The 1,4- hydroxy groups on the same 
side lead to much lower chemical shifts, 64.3 ppm for 
the Ch atom closer to the ribbon edge and 63.0 ppm for 
the other Ch atom. If wider 10-AGNR is used, the cor- 
responding chemical shifts are 67.1 and 64.5 ppm. For 
1,4-hydroxy groups on two opposite sides, the averaged 
chemical shifts are 63.4 and 65.4 ppm for 9-AGNR and 

10- AGNR, respectively. 

Therefore, isolated hydroxy group and hydroxy pair 
in many cases give Ch chemical shifts around 70 ppm. 
But in some other cases (1,4-hydroxyl), it is also possible 
to obtain chemical shifts between 60-70 ppm. At the 
fully-oxidized limit, hydroxyl groups will bond to two 
opposite sides for each neighboring carbon atom pair. 
However, there is still a large repulsion between hydroxy 
groups. The optimized C-C bond length is 1.64 A, 
which is much larger than a normal C-C bond length and 
leads to an energetically very unfavorable structure. The 
calculated Ch chemical shift is 91.7 ppm. 



C. Mixed epoxy and hydroxy groups 

In real GO samples, both epoxy and hydroxy groups 
exist, and they are expected to be randomly distributed 
on the graphene basal planes. Therefore, a realistic GO 
model must include both epoxy and hydroxy groups and 
consider the interaction between them. For simplicity, 
we still start from the lowest concentration limit using 
ribbon models. Obtained structure motifs can be used 
as building blocks for a more realistic GO model with 
higher oxidation group concentration. 

Based on energetics consideration, two stable local 
structures consisting both epoxide and hydroxyl have 
been proposed. [J^j Here we adopt them in our ribbon 
model and name them structure A and structure B (Fig- 
ure[3]), respectively. When structure A is put on 9-AGNR 
at a D-site (Figure [3^), the chemical shifts of two C^; 
atoms are close, with an average of 56.4 ppm. The dis- 
tance between the two epoxide carbon atoms is 1.47 A, 
slightly shorter than those of isolated epoxide. Chemical 
shifts for the two Ch atoms are 66.4 and 70.5 ppm, with 
the one closer to ribbon edge 4.1 ppm larger. This trend 
is in agreement with isolated 1,4-hydroxyl pair. When 
structure A is located at a P-site (Figure ISb) , the aver- 
age Ce chemical shift is 58.2 ppm. The optimized C-C 
bond length in the epoxide is still 1.47 A. The chemi- 
cal shifts for the two Ch atoms are 69.3 and 65.5 ppm, 
with the up Ch about 4 ppm larger. This is possibly due 
to the hydrogen bond formed between the two hydroxy 
groups. Similar results have been obtained for 10-AGNR. 
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For structure B at a D-site (Figure |3]:), chemical shifts 
of two Ce atoms are 62.2 and 69.1 ppm. While the chem- 
ical shifts of two Ch atoms are close, with an average of 
70.5 ppm. Here, the bond length between two epoxide 
carbon atoms is 1.49 A, and it is 1.56 A between the 
two hydroxyl carbon atoms. For P-site (Figure [3ji) , the 
chemical shifts of the two Ce and two Ch atoms are 
close, and their averages are 62.8 and 69.2 ppm, respec- 
tively. The distance between two epoxide carbon atoms 
is 1.48 A. Such a bond length slightly shorter than that 
of an isolated epoxide may be an indicator of lower chem- 
ical shift. In the 10-AGNR case, we get similar result for 
P-site, while for D-site, the chemical shift difference for 
two Ce atoms becomes much smaller, with an average 
chemical shift of 64.7 ppm. This result indicates that the 
big difference of the two C e chemical shifts observed in 
9-AGNR is not intrinsic. 



FIG. 3: Structure A locates at (a) a D-site or (b) a P-site in 
9-AGNR. Structure B locates at (c) a D-site or (d) a P-site 
in 9-AGNR. Two hydroxyl pairs in (e) zigzag or (f) armchair 
directions. Only part of the supercell used in calculations is 
shown. Carbon is in gray, oxygen in red, and hydrogen in 
white. 

Now, we have seen that, by adding proximate hydroxy 
groups, it is possible to decrease the chemical shift of 
epoxy group from about 70 ppm for isolated group to- 
wards the experimental value at 60 ppm. Therefore, in 
a realistic GO model, epoxy groups should be close to 
hydroxy and other epoxy groups. With higher oxidation 
group concentration, many GO structure models have 
been reported in the literature. [H, [12, [H, [l^ We choose 
some typical structures for NMR simulations. 

Based on a 2x2 supercell, Boukhvalov et al. [l^] pro- 
posed the most stable structure of GO with a 75% cov- 
erage of oxidation groups. Their model (Figure |4^) has 
a chemical formulism of C8(OH)40, and its unit cell is 
composed of a structure B unit and a 1,2-hydroxy pair. 
An important feature of the Boukhvalov model is that 
the hydroxy groups form chains, with hydrogen bonds 
forming within the chains, which is a main reason why 
such a structure is energetically very stable. The aver- 
age chemical shifts of C_e and Ch are 66.8 and 76.8 ppm, 
respectively. Both are much larger than their experimen- 
tal values. The Ce-Ce bond length is 1.51 A, and the 
averaged length of C-C bonds in hydroxy chains is 1.58 
A, larger than that of isolated hydroxy pair. We find 



that the hydroxyl molecules are the main reason of the 
wrinkling of carbon skeleton, as also suggested in a pre- 
vious study, dl] When hydroxy groups form chains, it 
elongates the corresponding C-C bonds, which may then 
generate larger Ch chemical shifts. 




FIG. 4: Typical GO structure model proposed by (a) 
Boukhvalov et al., (b) Lahaye et al., (c) Yan et al., and 
(d) Wang et al.. The supercells used in our calculations are 
marked with red lines. Carbon is in gray, oxygen in red, and 
hydrogen in white. 

With the same 2x2 supercell, calculations by Lahaye 
et al. suggest that too many hydroxy groups will lead to a 
too large tension. [2^ Therefore, instead of C8(OH)40, 
a GO model (Figure |3Jd) with a chemical composition 
C8(OH)202 was suggested. In this model, hydroxyl 
molecules are attached to carbon atoms directly adja- 
cent to epoxides, but at an opposite side of the carbon 
plane. Hydroxy groups do not form chains any more. 
Actually, this structure can be considered as a composi- 
tion of a structure A unit and an additional epoxy group 
per unit cell. After optimization, the length of the bond 
between two epoxide carbon atoms is 1.49 A, and the 
average length of the three C-C bonds connecting to a 
hydroxy group is about 1.54 A. Our NMR calculations 
give an average C^; chemical shift of 60.6 ppm, and the 
chemical shifts of the two Ch carbon atoms are both 
75.6 ppm. Compared to experiments, the 1,4-hydroxyl 
configuration adopted in this model leads to a too large 
chemical shift. 

Based on a systematic energetics analysis, Yan et al. 
argued that all structures with a negative formation en- 
ergy are fully oxidized GO with 1,2-hydroxyl pairs form- 
ing a chain-like structure and epoxy groups on remaining 
C atoms. 25] One example is C6(OH)202 (Figure Ht) , 
where 1,2-hydroxyl chains are along the armchair direc- 
tion instead of the zigzag direction as in the Boukhvalov 
model. The distance between two C^ atoms is about 1.49 
A, and the average Ch-Ch bond length is 1.58 A, similar 
to those in chains in the zigzag direction. The average 
Ce chemical shift is 61.9 ppm, close to the experimental 
value. While the average Ch chemical shift is still too 
large (78.5 ppm), since hydroxy groups form chains. 
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Another GO structure model is constructed based on 
the structure A building block, jl^] As shown in Figure 
mi, in this model, structure A units are connected along 
the armchair direction. Here, the ratio of C(sp2)/C(-0- 
)/C(-OH) is 1:1:1. The averaged C-C distance for epoxy 
groups is 1.49 A, while that for hydroxy pairs is 1.57 A. 
The average chemical shifts for C^; and Gh atoms are 
61.1 and 73.8 ppm, respectively. This structure model 
gives the best agreement with NMR experiments among 
the several models we studied here. We note several im- 
portant structure features of this structure: the proxim- 
ity between hydroxyl and epoxide, 1,2-hydroxyl pairs but 
without forming a chain structure, and a balanced ratio 
between hydroxyl and epoxide groups. 



D. sp carbon 



The chemical shift of sp^ carbon is also very important, 
since it is an indicator of the distribution of this carbon 
species. The chemical shifts of sp^ carbon in nanorib- 
bon models, including clean AGNRs, are generally be- 
tween 119-130 ppm, with most of them around 120-125 
ppm. They are smaller than the experiment value (129- 
133 ppm), which may be due to the limited width of 
the nanoribbons considered in this study and the hy- 
drogen atoms at the ribbon edges. For sp'^ carbon in 
an isolated C-C double bond surrounding by epoxy and 
hydroxy groups, we generally get a chemical shift much 
higher than the experimental value. All those GO models 
in Figure m which is not fully oxidized, have isolated C-C 
double bonds. The chemical shifts of isolated C-C double 
bonds are 143.4, 147.2, and 137.2 ppm for the model in 
Figure S^, |3)d, and |3Jl, respectively. Therefore, isolated 
sp^ carbon pair is not the main form of sp^ carbon in 
GO. 

Yan et al. have suggested that a stable GO structure 
is composed of fully-oxidized regions like that shown in 
Figure |3)d and sp^ carbon strips between them. 25| Al- 
though their fully oxidized GO model gives too high Ch 
chemical shifts, our calculated results support their de- 
scription on sp^ C. For aromatic strips between fully ox- 
idized regions, we obtain an averaged sp^ chemical shift 
about 134 ppm, [4l| which is in good agreement with ex- 
periments. We note that the large radius limit of the ^^C 
chemical shift of CNTs are much lower than 130 ppm. 
[Slj Therefore, very large graphene area is also not ex- 
pected to be existed in GO. This is consistent with the ex- 
perimental cross peak between sp^ C and epoxy/hydroxy 
groups in 2D NMR. [g^l Finally, we reach the following 
picture about sp^ carbon in GO: the main form of sp^ 
carbon is small size aromatic clusters between highly ox- 
idized regions. 



E. Other oxidation groups 

The assignment of the three small peaks in GO NMR 
spectrum is more difficult than that of the three strong 
resonances. Only recently, based on ^H-^'^C cross po- 
larization (CP) spectrum, it was suggested that the 101 
ppm signal is resulted from non-protonatcd carbons, [2^] 
possibly a peripheral structure of GO containing five- 
and six-membered-ring lactols. Q The 167 and 190 ppm 
peaks have been tentatively assigned to ketone and ester 
carbonyl, respectively. '9, 2ll] 







FIG. 5: (a) An AGNR model with 5- and 6-membered-ring 
lactols, carboxylic acid groups, and ketone groups. AGNR 
models with (b) an epoxy pair and (c) an epoxy-hydroxy 
pair. Two dimensional models with (d) epoxy pair and (e) 
epoxy-hydroxy pair. The unit cells used in our calculations 
are marked by the red line. Carbon is in gray, oxygen in red, 
and hydrogen in white. 

To confirm those assignments reported in the litera- 
ture, we construct an AGNR model (Figure [5^), which 
contains 5- and 6-membered-ring lactols, esters, car- 
boxylic acid, and ketone groups. Chemical shifts of the 
sp'^ carbon atom in 5- and 6-membered-ring lactols are 
104.9 and 96.9 ppm, respectively. Both are close to 101 
ppm. The averaged chemical shift of C=0 is 191.6 ppm, 
which is also in agreement with the 190 ppm peak in pre- 
vious assignments. However, we note that the chemical 
shift of C=0 can be strongly affected by its environment. 
For example, in Fig[5t, the edge C=0 group has a chem- 
ical shift only 173.1 ppm. The chemical shift of the edge 
carboxyl carbon is 162.2 ppm, close to 167 ppm peak in 
experiments. There is also a COOH structure in 5- and 
6-membered-ring lactols, which gives chemical shifts of 
166.3 and 163.2 ppm respectively. Therefore, our calcu- 
lations generally confirms previous assignment. 

However, in a previous 2D NMR experiment, f23| cross 
peak has been observed for the 101 ppm peak, but not the 
other two minor peaks. Therefore, the 101 ppm peak is 
different in the three minor peaks. It may comes from on- 
plane groups instead of edge group. In a previous study, 
we have suggested that epoxy pair and epoxy-hydroxy 
pair may exist in highly oxidized GO samples. i26t] It is 



6 



interesting to see what is their NMR signals. We use both 
ribbon model and 2D model to simulate their chemical 
shifts. Figure [5b shows an AGNR with an epoxy pair, 
where both carbon atoms bound by the epoxy pair have 
a chemical shift around 103.7 ppm. The chemical shift 
for epoxy-hydroxy pair in the ribbon model (Figure [St) 
is 105.9 ppm. Bulk models give similar results, with a 
101.4 ppm average chemical shift for epoxy pair (Figure 
[5]i) and 106.8 ppm for epoxy-hydroxy pair (Figure [5]3) . 
Therefore, the NMR signal at 101 ppm strongly support 
the existence of epoxy pair we proposed in a previous 
cutting mechanism study [13] and XPS simulation. [2^ 
Although its chemical shift is also close to 101 ppm, our 
results can not be used to support the existence of the 
epoxy-hydroxy pair, since the 101 ppm signal disappears 
in CP spectrum. 

F. Discussion 

Isolated epoxide presents a Ce chemical shift around 
70 ppm. It decreases to around 60 ppm in most cases, 
when epoxide is in close proximity with hydroxyl. A pre- 
vious study suggested that the sole presence of 1,2- 
ether epoxy group is not stable due to the accordingly 
created big tension on the carbon grid, and hydroxyl is 
needed to safe guard the stability of the structure. There- 
fore, the proximity of epoxy and hydroxy groups has both 
significant geometric effects and electronic effects. It's in- 
teresting to see whether the change of the chemical shift 
of the epoxy group is mainly a geometric or an electronic 
effect. 

For this purpose, a computer experiment has been per- 
formed by calculating chemical shifts of an artificial sys- 
tem, which has the hydroxy groups taken off from a previ- 
ously optimized geometry but fixing all other atoms. For 
the structure shown in Figure the average C^; chem- 
ical shift becomes 76.9 ppm after hydroxyl groups being 
taken off. For the structure in Figure l^b, the averaged 
Ce chemical shift also becomes 68.8 ppm, with obvious 
increase. Therefore, electronic effect plays an important 
role in the hydroxy induced chemical shift decrease of 
epoxide. 

Our results suggest that hydroxy chain is not widely 
existed in GO, since its chemical shift is too high. To 
confirm this conclusion, we further consider the effect 
of nearby epoxy groups. Similar calculations have been 
performed by removing epoxy groups from the optimized 
geometries with hydroxy chains. For the structure shown 



in Figure [4^, the average Ch chemical shift is 78.2 ppm 
after removal of epoxide, which is only 1.4 ppm larger 
than the original value. For the structure in Figure 011, 
we also only get a 1.7 ppm chemical shift difference by 
removing the epoxy groups. 

Therefore, the great geometrical distortions brought 
by the hydroxy chains should be the main reason of the 
chemical shift increase related to isolate hydroxy pairs. 
To see if there is any other effect, we consider two neigh- 
boring hydroxy pair. They can arrange along the zigzag 
direction (Figure [3^) or along the armchair direction 
(Figure [3]:) . In both cases, the Cr-Gh bond length is 
about 1.55 A, much smaller than those in hydroxy chains. 
Therefore, the geometrical effect is small in these two 
cases. 

In the zigzag direction case, the chemical shifts are 
74.5, 67.2, 72.9, and 71.6 ppm for Gh carbon atoms from 
left to right, respectively. In the armchair case, the corre- 
sponding chemical shift is 74.1, 70.1, 70.4, and 74.2 ppm, 
from up to down. We can clearly see that the hydroxy 
group with a hydrogen bond pointed to it has a higher 
Ch chemical shift. Therefore, hydrogen bonding may 
also be an important factor in increasing the chemical 
shift of hydroxy group in hydroxy-pair chain. 



IV. CONCLUSIONS 

Our results are consistent with the Lerf model of GO, 
where the main species in GO are epoxide, hydroxyl, 
and sp^ carbon. Their chemical shifts are very sensitive 
to their chemical environment, which leads to the result 
that peaks in GO NMR spectrum are very broad. When 
it is isolated, the chemical shift of epoxide is too high. 
Therefore, epoxide prefers to in close proximity with hy- 
droxyl. Chemical shift of hydroxyl groups is closely re- 
lated to its geometrical structure. Isolated 1,2-hydroxy 
pairs have stable chemical shifts around 70 ppm. How- 
ever, when they form hydroxy chains, the chemical shift 
becomes much higher mainly due to the elongated C- 
C bonds. As a result, hydroxy chains widely proposed 
in literature based on energetics considerations should 
not be an important structure motif of GO. Our calcu- 
lations also confirmed the existence of aromatic carbon 
clusters among highly oxidized regions. Considering its 
cross peak in 2D NMR, the small peak at 101 ppm is more 
likely contributed by on-plane groups, such as epoxy pair 
we proposed earlier. 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 

Zhang, S. V. Dubonos, I. V. Grigorieva and A. A. Firsov, 

Science 306, 666 (2004). 
[2] A. K. Geim and K. S. Novoselov, Nature Mater. 6, 183 

(2007). 

[3] M. I. Katsnelson, Mater. Today 10, 20 (2007). 



[4] A. H. Castr Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

[5] D. R. Dreyer, S. Park, C. W. Bielawski and R. S. Ruoff, 

Chem. Soc. Rev. 39, 228 (2010). 
[6] S. Stankovich, D. A. Dikin, G. H. B. Dommett, K. M. 



7 



Kohlhaas, E. J. Zimney, E. A. Stach, R. D. Finer, S. T. 

Nguyen and R. S. RuofT, Nature 442, 282 (2006). 
[7] S. Gilje, S. Han, M. Wang, K. L. Wang and R. B. Kaner, 

Nano Lett 7, 3394 (2007). 
[8] S. Park and R. S. Ruoff Nature Nanotechnology 4, 217 

(2009). 

[9] W. Gao, L. B. Alemany, L. Ci and P. M. Ajayan, Nature 

Chemistry 1, 403 (2009). 
[10] D. A. Dikin, S. Stankovich, E. J. Zimney, R. D. Piner, 

G. H. B. Dommett, G. Evmenenko, S. T. Nguyen and R. 

S. Ruoff, Nature 448, 457 (2007). 
[11] X. Wu, M. Sprinkle, X. Li, F. Ming, C. Berger and W. 

A. de Heer. Phys. Rev. Lett. 101, 026801 (2008). 
[12] L. Wang, K. Lee, Y.-Y. Sun, M. Lucking, Z. Chen, J. J. 

Zhao and S. B. Zhang, ACS Nano, 3, 2995 (2009). 
[13] G. M. Scheuermann, L. Rumi, P. Steurer, W. Bannwarth 

and R. Miilhaupt, J. Am. Chem. See. 131, 8262 (2009). 
[14] G. Eda, Y.-Y. Lin, C. Mattevi, H. Yamaguchi, H.-A. 

Chen, L-S. Chen, C.-W. Chen and M. ChhowaUa, Adv. 

Mater. 22, 505 (2010). 
[15] I. V. Lightcap, T. H. Kosel and P. V. Kamat, Nano Lett 

10, 577 (2010). 

[16] M. Mermoux, Y. Chabre and A. Rousseau, Carbon 29, 
469 (1991). 

[17] H. He, T. Riedl, A. Lerf and J. Khnowski, J. Phys. Chem. 

100, 19954 (1996). 
[18] A. Lerf, H. He, T. Riedl, M. Forster and J. Khnowski, 

Sohd State Ionics 101-103, 857 (1997). 
[19] A. Lerf, H. He, M. Forster and J. Khnowski, J. Phys. 

Chem. B 102, 4477 (1998). 
[20] W. Cai, R. D. Piner, F. J. Stadermann, S. Park, M. A. 

Shaibat, Y. Ishii, D. Yang, A. Velamakanni, S. J. An, M. 

StoUer, J. An, D. Chen and R. S. Ruoff, Science 321, 

1815 (2008). 

[21] T. Szabo, O. Berkesi, P. Forgo, K. Josepovits, Y. Sanakis, 

D. Petridis and I. Dekany, Chem. Mater. 18 2740 (2006). 
[22] D. W. Boukhvalov and M, L Katsnelson, J. Am. Chem. 

Soc. 130, 10697 (2008). 
[23] R. J. W. E. Lahaye, H. K. Jeong, C. Y. Park and Y. H. 

Lee, Phys. Rev. B 79 125435 (2009). 
[24] Z. Li, W. Zhang, Y. Luo, J. Yang and J. G. Hon, J. Am. 

Chem. Soc. 131, 6320 (2009). 



[25] J. A. Yan, L. Xian and M. Y. Chou, Phys. Rev. Lett. 

103, 086802 (2009). 
[26] W. Zhang, V. Carravetta, Z. Li, Y. Luo and J. Yang, J. 

Chem. Phys. 131, 244505 (2009). 
[27] X. Gao, L. Wang, Y. Ohtsuka, D. Jiang, Y. Zhao, S. Na- 

gase and Z. Chen, J. Am. Chem. Soc. 131, 9663 (2009). 
[28] K. N. Kudin, B. Ozbas, H. C. Schniepp, R. K. 

Prud'homme, L A. Aksay and R. Car, Nano Lett. 8, 36 

(2008) 

[29] M. A. L. Marques, M. d'Avezac and F. Mauri, Phys. Rev. 

B 73, 125433 (2006). 
[30] E. Zurek, C. J. Pickard, B. Walczak and J. Autschbach, 

J. Phys. Chem. A 110, 11995 (2006). 
[31] L. Lai, J. Lu, W. Song, M. Ni, Lu Wang, G. Luo, J. Zhou, 

W. N. Mei, Z. Gao and D. Yu, J. Phys. Chem. C 112, 

16417 (2008). 

[32] E. Zurek, C. J. Pickard and J. Autschbach, J. Am. Chem. 

Soc. 129, 4430 (2007). 
[33] E. Zurek, C. J. Pickard and J. Autschbach, J. Phys. 

Chem. C 112, 11744 (2008). 
[34] E. Zurek, C. J. Pickard and J. Autschbach, J. Phys. 

Chem. A 113, 4117 (2009). 
[35] M. D. Segah, P. J. D. Lindan, M. J. Probert, C. J. 

Pickard, S. J. Hasnip and M. C. Payne, J. Phys.: Con- 
dens. Matter 14, 2717 (2002). 
[36] S. J. Clark, M. D. Segall, C. J. Pickard, P. J. Hasnip, M. 

J. Probert, K. Refson and M. C. Payne, Z. Kristallogr. 

220, 567 (2005). 
[37] D. Vanderbih, Phys. Rev. B 41, 7892 (1990). 
[38] Y. Zhang and W. Yang, Phys. Rev. Lett. 80,890 (1998). 
[39] C. J. Pickard and F. Mauri, Phys. Rev. B 63, 245101 

(2001). 

[40] J. R. Yates, C. J. Pickard and F. Mauri, Phys. Rev. B 

76, 024401 (2007). 
[41] See EPAPS supplementary material at 
|http://dx. doi.org/10.1063/ for optimized geometries and 

detailed chemical shift values. 
[42] Y. W. Son et al., Phys. Rev. Lett. 97, 216803 (2006) 
[43] L. Nu, Z. Li, and J. Yang, J. Phys. Chem. C 113, 16741 

(2009) 



